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Abstract 



We present a microscopic theory for the low temperature metamagnetic 
phase diagram of H0M2B2C that agrees well with experiments. For the same 
model we determined the zero field ground state as a function of tempera- 
ture and find the c-axis commensurate to incommensurate transition in the 
expected temperature range. The complex behaviour of the system originates 
from the competition between the crystalline electric field and the Ruderman- 
Kittel-Kasuya-Yosida interaction, whose effective form is obtained. No es- 
sential influence of superconductivity has to be invoked to understand the 
magnetic phase diagram of this material. 
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The recent interest in H0M2B2C and similar borocarbide compounds is motivated by 
the possibility of a detailed study of the mutual interaction between superconductivity (SC) 
|]J and magnetic order (MO) |2|,[J coexisting in a few of these materials as bulk proper- 
ties. The superconducting critical temperature for HoNi 2 B 2 C is T c = 8 K [0 and the upper 
critical field is about 2 ~ 3.5 kG ||. Temperature dependent measurements show a 
pronounced anomaly of H^(T) around 5 K which nearly leads the material to reentrance 
into the normal state fHH. In the same temperature range several magnetically ordered 
structures are observed in the Ho f-electrons sub-lattice: a commensurate (C) antiferromag- 
netic phase below 5 K @,[|, an incommensurate (IC) c-axis complex spiral state (5 K < T 
< 6 K) [0,||[7| and an a-axis IC modulation in a narrow range of temperature around 5.5 
K Although no satisfactory theory is available presently, many experiments point to 

a correlation between this complex magnetic phase diagram and the SC anomalies p|,|T0 



Much insight on the magnetism of this material can be gained from a number of experiments 



reporting metamagnetic transitions at low temperature and fields higher than P,|TT||. In 
particular the detailed anisotropic metamagnetic phase diagram at T = 2 K presented in 
Ref. [[E|, can be used to extract many features of the magnetic interaction. In this work 
we propose a realistic microscopic model for the Ho 4 /-electrons sub-system of H0M2B2C 
which reproduces the main features of the low temperature metamagnetic phase diagram as 
well as the zero field sequence of phases as function of temperature. It is important to note 
that many physical elements need to be included in the present theoretical description to 
have reasonable agreement with experiments and likely a model of a comparable complexity 
is needed to treat the mutual interaction between SC and MO. 

The coexistence and weak coupling of the two phenomena of SC and MO is due to the 
different degrees of localisation of electrons in the borocarbides. LDA calculations show 



that the conduction band is composed mainly of Ni 3<i-electrons ||13|| , which undergo the 
superconducting transition. On the other hand magnetic properties are related to the well 
localised electrons in the incomplete 4 /-shell of Ho. The exchange interaction between these 
two electron systems is mediated by the small fraction of Ho 6s and 5d character in the 



conduction band. As far as magnetic properties are concerned the conduction electrons can 
be eliminated in a standard way leading to an effective RKKY exchange interaction among 
the stable Ho 4/ moments. The appropriate Hamiltonian is then given by [|TJ]: 

H = E[K/(Ji) - to ■ B] - ~ £ J{i, j) J, • Jj (1) 

i ij 

This Hamiltonian include the crystal electric field (CEF) single ion part 7Y c /(Jj) expressed in 
terms of the total angular momentum Jj, the Zeeman interaction between the local magnetic 
induction B and the magnetic moment /Xj = /^BgJi (//b is the Bohr magneton and g — | 
the gyromagnetic ratio for Ho), and the effective RKKY exchange interaction. The direct 



dipole-dipole interaction is not relevant in borocarbides [15] and this allows us to use a 



magnetic induction field B independent of the position. The CEF single ion Hamiltonian 



we use is the one extracted from neutron diffraction experiments in Ref. [|16j and contains 
no adjustable parameter. Its ground state is a T 4 singlet and the first exited states are a 
doublet at 0.15 meV from the ground state and a I\ singlet at 0.32 meV. The other 13 CEF 
states have much higher excitation energies (> 10 meV) and their matrix elements with 
the low energy quartet is very small, therefore they may be neglected in the whole range of 
temperature and fields which we explored. Regarding the magnetic interaction no previous 
knowledge is available about the RKKY function J~(i,j) and an important aim of this work 
is to obtain a realistic model for it. 

Since we are dealing with a three dimensional system of large angular momenta (Jho = 
8) it is possible to treat its Hamiltonian at a mean field (MF) level. Introducing the mean 
thermal value (Jj) and neglecting the terms containing two sites fluctuations it is possible to 
decouple the dynamics of the different sites. The single ion MF Hamiltonian is then given 
by 0: 

H mf (i)=H cf (J i )-J i -B e i +£ (2) 

in which B^ = (/i^B + J2j<?(hj) (Jj)) is the effective molecular field and £ = 
\ {^i)Yjj <J{h j) (Jj)- The diagonalisation of this single site Hamiltonian can be easily 



achieved numerically. The single ion Gibbs free energy density Ff^ and the corresponding 
average angular momentum ( Jj) can be computed as functions of B ? ? and T. This leads to the 
self consistent MF equations for the (Jj) which can be solved iteratively. Two experimental 
evidences can be used to reduce the number of independent sites whose MF magnetisation 
should be computed. First of all, almost all the observed structures in this compound share 
the property of ferromagnetic alignment in the ab plane. The only exception is the a-axis 
modulation whose structure is not yet clear. However it appears not to coexist microscop- 
ically with the c-axis ones and we will neglect it. Therefore we impose ferromagnetic 
alignment in the plane from the beginning thus reducing the calculation to one dimension. 
Then the RKKY interaction have only c-dependence and, taking the magnetic moment in 
the site as the reference, it may be parametrised with the help of Ji = Y^ua <?(®,ji), with 
jj running on all the sites of the i th plane. The Ji with % > 1 are the interaction of the 
reference moment with those in the neighbouring layers and Jo is the interaction with the 
other moments in the same plane. The second simplification is related to the CEF structure, 
in fact the c-axis is a very hard direction in the whole range of temperature we are interested 
in, therefore we force the moments to lie in the ab plane neglecting their small out-of-plane 
component. 

In order to establish the actual MF ground configuration for the (Jj) out of the possible 
stable ones we introduce the Helmholtz free energy (HFE) density given by: 

F(r) = F M (B) + f - ^ = F M (B) - 2,M 2 - !f (3) 



with B(r) = H(r) + 47rM(r) JT7J. The integral of F(r) over all space is the proper ther- 
modynamic function to be minimised when external fields are kept constant ||17|| . The 
magnetisation contribution (~ 10 kG) is not negligible with respect to the typical external 
field (4-25 kG). However the full problem of solving for M(r) in a finite sample with a given 
geometry is beyond our purpose. This is indeed a typical problem in thermodynamics of 



magnetic materials jI7] whose most practical solution is to assume cylindrical symmetry 
around the external field H in order to eliminate the spatial dependence of both B and M. 



The homogeneous magnetisation of the sample is calculated as M = ^-(Ji) where V c = 
65 A 3 is the volume of the unit cell and the bar indicates the average on all ions. Similarly 
the contribution of the magnetisation to the HFE per unit cell can be written as 2-kM 2 V = 

Until now the actual RKKY interaction among magnetic moments remained unspecified. 
In order to achieve a convenient parametrisation for it we make an extensive use of the 



T = 2K magnetisation data in Ref. |12| . From the clear presence of flat magnetisation 
plateaux and from the value of the magnetisation as a function of the angle we argue that 
the magnetic moments of the ions are almost at the saturation value J s = J# = 8 and 



they are locked in one of the four equivalent (110) in-plane easy directions | 12| . The sharp 
metamagnetic transitions are then due to first order transitions among different arrangement 
of the moments along the c-axis. We assume that the relevant phases in the field-angle phase 
diagram are the easiest commensurate structures with the observed magnetisation (they are 
listed in the upper part of table |) . We notice that all the phases in table | may be represented 
in a chain of six unit cells with periodic boundary conditions. This is important because the 
Fourier transform of the interaction for a system with six sites have only four free parameter, 
namely Ji with i = 0, 1, 2, 3. Any further Fourier component may be removed with a 
renormalisation of these parameters (i.e. the is equivalent to Jq). Neglecting the CEF 
energy and entropy it is possible to calculate analytically the HFE for each configuration. 
Some of these energies are listed in table |. From the two 9 = 0° metamagnetic transitions 
it is possible to establish two conditions for the parametrisation of the RKKY interaction. 
Following the convention of Canfield et. al. [0 we call the metamagnetic transition fields 
H c i = 4.1 kG and H C 2 = 11.1 kG (this value is somewhat higher than the 10.6 kG given in 



Ref. fl2[]), they should not be confused with the superconducting critical fields. Imposing 
the energy of the AF2 and AF3 phases to be equal at H cl and the ones of AF3 and P at H c2 
we obtain the following relations: 

Jx = -J* ~ ^fH c2 - \e m (4) 



Jz = J2 + ^ry-{H C 2 — H c \) + -E M 
bJ s o 



where E M = 2 " ( f a) = 8.1 lO" 3 meV comes from the contribution of the magnetisation 
to the HFE. They assure that the relative stability of the three phases for 9 = 0° is the 
observed one. In order to use the easiest possible model we set J$ = 0, which implies J\ = 
—8.0 10~ 3 meV and J2 = —2.4 10™ 3 meV. The final freedom in the model is the parameter 
Jo which works basically as a self interaction and cannot be extracted starting from the 
energy differences among magnetic structures. We use J§ = 4.8 10~ 3 meV in order to have 
the transition between the paramagnetic and the incommensurate state at the experimental 
value T = 6 K ||. Moreover this choice make the internal molecular field large enough (~ 
15 kG) to maintain the moments close to the saturation regime, as required. 

After the model has been defined and all its parameters fixed we present now the numer- 
ical results of the complete self consistent MF calculation for the field-angle phase diagram 
at T = 2 K. The starting values for (Jj) in the iteration algorithm are a set of random num- 
bers and the possible periods allowed are in the range from one to nineteen planes. Typical 
magnetisation curves are evaluated and they are shown in fig. [I], they refer to experimental 
geometries with different angles 9 between the external field H and the closest magnetic easy 
axis (110). The resulting metamagnetic phase diagram is presented in fig. |2|. It contains 
all the phases in the table | and it agrees remarkably well with the experimental one. In 
particular all the transition lines show simple trigonometrical dependence as a function of 
the angle 9 as well as the corresponding magnetisation values at the plateaux. In all the 
regions where the stable phases are AF2, AF3, F3 and P we observe remarkable quantita- 



tive agreement with the experimental data in Ref. | 12| . The main qualitative differences 
between our model and experiments is the presence of two additional phases, the F2 and 
the C6, whose magnetisation is different from the reported ones. However the strongest dis- 
agreement is in the region of relatively low field and high angles where experimentally the 
strongest hysteresis is found. Another minor point is the absence of direct AF3-P transition 
for 9 7^ 0°. Experimentally the direct transition is observed for small angles up to 9 = 6° in 
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the Ho compounds, but seems to have a much wider range in the similar phase diagram of 
DyNi 2 B 2 C ||18|| . We would like to stress however that it is not possible to improve this phase 
diagram simply refining the parameters for the RKKY interaction. For example allowing 
for J~3 ^ it is possible to stabilise at low field an additional phase, which is a distorted 
helix structure with wave length five, but no change appears in the phase transitions among 
the other phases. In addition, as explained before, the effect of further effective interactions 
with neighbouring planes beyond the third one may be eliminated by proper renormalisa- 
tion of the J7i with i < 3. Therefore if the two phases F2 and C6 are not observed in the 
experiments we have to conclude that additional interactions (i.e. magnetoelastic couplings) 
have to play a role in the stability of magnetic phases in HoNi 2 B 2 C. 

Starting from the same model it is possible to analyse the zero field behaviour as a 
function of temperature. In principle this requires some attention since the magnetic system 
is now embedded in a superconducting material. This implies important changes in the 



q ~ region of the Fourier transform of the RKKY function but leaves the relevant 
q ~ 7r region almost unchanged. Relying on this fact and on experiments on doped non- 
superconducting materials such as HoNi 2 _ x Co x B 2 C which show a magnetic behaviour very 



similar to the undoped superconducting one POinj , we will use our purely magnetic model 



for the description of the zero field phase diagram. At MF level the second order phase 
transition between the paramagnetic state and an ordered structure is expected to occur a 
the Q vector for which the J7"(q) has its maximum. In our model this correspond to Q c = 
0.78 7r in the (001) direction, not far from the experimental value Q exp = 0.9l7r 0. The 
helical state is preferred with respect to the longitudinal modulated structure due to the ab 
easy plane for the moments. This truly incommensurate structure can be the ground state 
of the system only as long as the average moment per ion is small enough, i.e. close to the 
transition temperature. Lowering the temperature the ordered state develops and the CEF 
part of the HFE, proportional to fourth and the sixth powers in (Jj), force the structure 
to find a commensurate compromise. Because in the self consistent MF treatment it is not 
possible to treat at the same time truly incommensurate structures and the CEF, we cannot 



observe the actual C-IC transition. However at a temperature of T = 5.5 K is observed 
a first order transition from AF2 to a helical state of wave length 17, where the moments 
point no longer only along the easy directions. This is the clear indication that the RKKY 
energy starts to become dominant with respect to the CEF potential and drives the system 
into a state whose wave number Q is closer to the maximum of the RKKY function. To 
obtain better quantitative agreement for the ordering wave number and for the temperature 
interval in which the incommensurate state is stable, the function J7"(q) has to be refined in 
the Q c region by including further parameters. 

In conclusion, we presented a microscopic model for the the rare earth borocarbide 
system HoNi 2 B 2 C which explains the main reported features of the anisotropic magnetic and 
temperature phase diagrams. The minimal model to achieve a semiquantitative description 
of the complex magnetic behaviour of the system needs to include realistic CEF and effective 
RKKY interaction among the planes. On the other hand no influence of SC needs to be 
included in the determination of the magnetic structures observed. 

We would like to thank W. Henggeler et al. for the data on the CEF states. A. A. would 
also like to thank M. Laad, P. De Los Rios and B. Canals. This work was performed under 
DFG Sonderforschungsbereich 463. 
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TABLE I. Stacking sequence of ferromagnetically ordered a6-planes along the c-axis for the 
phases found in the T = 2 K magnetic anisotropic phase diagram in fig. 2. At low temperature 
the CEF forces the moments to lie in one of the four easy direction (110) indicated by arrows. The 
external field forms an angle with the (-1,1,0) (\) direction. The third column is the relevant 
part of the HFE per site, the term — J] & — |^ may be added to have the total HFE. Hy = g^sH 
cos(#) is the projection of the field along the easy axis and all other symbols are explained in the 
text. We only give the HFE for the phases needed to compute the relations (4). 
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FIG. 1. Magnetisation vs. magnetic field for some representative angles 9 between the field 
and the closest (110) direction. For 6=0° the thick dotted line represent experimental data taken 
from Canfield et al [12]. and the thick continuous line give results of our calculation. The two 
experimental parameters entering our model are the transition fields of the two metamagnetic 
transition H c i = 4.1 kG and H C 2 = 11.1 kG. This figure should be compared with fig. 1 (a) of [12] 
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FIG. 2. Field vs. angle metamagnetic phase diagram for T = 2 K. The filled dots are calculated 
within our model and indicate the phase transitions where a sensible jump in the magnetisation is 
seen (AM > .2 ^b/Ho). The thick dashed lines are given in [12] as the best fit of the experimental 
data. Their functional form is H° 3 /sin(45° - 9), H° 2 /cos(45° - 6) and B° cl /cos6. The values we used 
for the proportionality constants are H° 3 = H° 2 = 8.4 kG and H° x =4.1 kG. The two additional 
thin curves refer to the phase boundaries not yet observed in the experiments. Empty symbols 
refer to the stability of the AF3 phase with respect to AF2 (circles) and F2 (diamonds). The 
apparent tetracritical point in the phase diagram is composed of two very close usual tricritical 
points (inset). 
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